Radial and elliptic flow in Pb+Pb collisions at the Large Hadron Collider 

from viscous hydrodynamics 
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A comprehensive viscous hydrodynamic fit of spectra and elliptic flow for charged hadrons and 
identified pious and protons from Au+Au collisions of all centralities measured at the Relativistic 
Heavy Ion Collider is performed and used as the basis for predicting the analogous observables for 
Pb+Pb collisions at the Large Hadron Collider at ^/s = 2.76 and 5.5 A TeV. Comparison with recent 
measurements of the elliptic flow of charged hadrons by the ALICE experiment shows that the 
model slightly over-predicts the data if the same (constant) specific shear viscosity 77/s is assumed 
at both collision energies. In spite of differences in our assumptions for the equation of state, 
the freeze-out temperature, the chemical composition at freeze-out, and the starting time for the 
hydrodynamic evolution, our results agree remarkably well with those of Luzum [M. Luzum, Phys. 
Rev. C 83, 044911 (2011)], indicating robustness of the hydrodynamic model extrapolations. Future 
measurements of the centrality and transverse momentum dependence of spectra and elliptic flow 
for identified hadrons predicted here will further test the model and shed light on possible variations 
of the quark-gluon transport coefficients between RHIC and LHC energies. 
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I. INTRODUCTION 



The first measurement of elliptic flow in Pb-|-Pb colli- 
sions at the Large Hadron Collider (LHC) has just been 
reported [1]. The elliptic flow coefficient V2 character- 
izes the momentum anisotropy of final particle emission 
in non-central heavy-ion collisions relative to the event- 
plane, defined by the beam direction and the minor axis 
of the nuclear overlap region in the collision. It describes 
the efficiency of the medium generated in the collision to 
generate from an initial spatial deformation of its density 
distribution an asymmetry in the final momentum distri- 
bution, through interactions of the medium constituents. 
This efficiency increases with the coupling strength be- 
tween those constituents and becomes maximal for an in- 
finitely strongly coupled medium. In this limit the mean 
free path of the constituents becomes as small as allowed 
by the uncertainty relation [2], and the medium reaches 
very quickly a state of approximate local thermal equi- 
librium which allows to describe its evolution with fluid 
dynamics. For given initial spatial deformation of the col- 
lision fireball, ideal fiuid dynamics (which assumes zero 
mean free path) is expected to generate the largest possi- 
ble elliptic flow [3] . Shear viscosity, a consequence of non- 
zero mean free paths and limited from below by quantum 
mechanics [2, 4], will lead to a suppression of V2 [5, 6]. 
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Compelling evidence for fiuid dynamical behavior of 
the collision fireballs created in ultrarelativistic heavy- 
ion collisions, with a very small ratio of shear viscosity 
to entropy density, ry/s, has been found in heavy-ion col- 
lisions at the Relativistic Heavy Ion Collider (RHIC) [7- 
10]. The new data from the LHC confirm this picture 
[1, 11, 12] and agree, at least qualitatively, with hydro- 
dynamic predictions of elliptic flow for Pb+Pb collisions 
at the LHC [13-19]. 

The purpose of the present article is to explore how 
good this agreement is quantitatively, and to what ex- 
tent the present and future LHC elliptic flow data can 
tell us novel facts about the transport behavior of hot 
QCD matter at temperatures that exceed those accessi- 
ble at RHIC but are within reach at the LHC. Similar to 
the analyses in [11, 14, 16, 19], but different from recent 
hybrid model studies in [18, 20], we base our analysis on 
a purely hydrodynamic approach. While this ignores the 
fact that the late dilute hadronic stage of the expansion 
is very dissipative and not well described by fluid dynam- 
ics (neither ideal [21] nor viscous [22]), the importance 
of the hadronic phase for the development of elliptic flow 
is expected to be much reduced at the LHC relative to 
RHIC [18, 23]. As in Refs. [11, 16, 19, 20], but differ- 
ent from Refs. [14, 18], we use viscous hydrodynamics 
with a non-zero (but constant, i.e. temperature inde- 
pendent) specific shear viscosity 77/s, adjusted to spectra 
and elliptic flow measurements at RHIC. Our fitted value 
77/5 = 0.20 (for CGC initial conditions, see below) is 25% 
larger than that used by Luzum and Romatschkc [11, 16] 
but agrees well with the value for the quark-gluon plasma 
(QGP) viscosity (77/s)qgp recently extracted from RHIC 
data by using a hybrid viscous hydrodynamic + Boltz- 
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mann approach (VISHNU [10, 22]). Calculations with such 
a hybrid approach arc numericaUy much more demanding 
than purely hydrodynamic simulations; a generalization 
of the present analysis using VISHNU will follow soon and 
should further improve the reliability of the LHC predic- 
tions. 



II. HYDRODYNAMIC FIT OF RHIC AU+AU 
DATA 

In this work, we employ (2-(-l)-d viscous hydrodynam- 
ics [24] with the lattice QCD based equation of state 
s95p-PCE [25, 26], which accounts for chemical freeze-out 
before thermal decoupling at Tchom = 165 MeV, to sim- 
ulate the expansion of the collision fireball. From the 
analysis [10] of charged hadron spectra and elliptic flow 
in 200 A GeV Au+Au collisions at RHIC we take over 
a value of 7^/5 = 0.20 (corresponding to MC-KLN initial 
conditions, see below) for the effective specific shear vis- 
cosity of the strongly interacting fluid. Using the insights 
obtained from the systematic parameter study presented 
in [26] , we initialize the hydrodynamic expansion at time 
tq = 0.6fm/c and decouple at Tdoc = 120MeV at both 
RHIC and LHC energies. For Au-|-Au collisions at RHIC 
energies these parameters allow for a good global descrip- 
tion of the hadron pj^-spectra and differential elliptic flow 
(see below). Lacking strong theoretical or phenomeno- 
logical guidance how to adjust their values for Pb-fPb 
collisions at the LHC, we here decided to keep them un- 
changed. 

At thermalization time tq, we assume that the 
shear stress tensor is given by its Navier-Stokes value 
tt'^'^ = 2r]a'^'^ (where a'^'^ =V^^u'^^ is the symmetric and 
traceless velocity shear tensor), and that the initial ex- 
pansion flow is entirely longitudinal with Bjorken profile 
Vz = z/t and zero transverse flow velocity. In Milne co- 
ordinates (r, X, y, if) this corresponds to an initial flow 
4- velocity {u^ , , m^, u^) = (1, 0, 0, 0). Kinetic freeze-out 
is implemented by converting the hydrodynamic output 
to particle spectra with the Cooper-Frye prescription [27] 
on a decoupling surface of constant temperature T^ec- We 
use a quadratic ansatz 5f{x,p)^p^p'^'!i^^{x) [6] for the 
viscous deviation from local thermal equilibrium of the 
local phase-space distribution function on the freeze-out 
surface. Our final hadron spectra include decay products 
from strong decays of all particles and resonances up to 
2 GeV mass [28] , using the resonance decay program from 
the AZHYDRO package.^ 

For the initial density profile we here use a Monte- 
Carlo version [29, 30] of the Kharzeev-Levin-Nardi [31] 
model (MC-KLN).^ The specific implementation used in 
this work is described in [18, 32]. The model gives for 



1 AZHYDRO is available at the URL 
http : //www. physics . ohlo-state . edu/"f roderma/. 

^ The Monte Carlo code is available at URL 



each event the gluon density distribution in the trans- 
verse plane. We assume it to thermalize by time tq, 
and convert the gluon density into entropy density [21]. 
Over one million Monte Carlo events are re-centered to 
the beam axis and rotated in the transverse plane such 
that their minor axis aligns with the impact parameter 
(i.e. their "participant plane" coincides with the reac- 
tion plane). After sorting them into collision central- 
ity bins according to their number TVpart of participat- 
ing ("wounded") nucleons, we average them to obtain 
a smooth average initial entropy density which is then 
converted to energy density using the equation of state. 
Using this smooth energy density as weight, we com- 
pute the initial eccentricity e — ^^^^^^^ and overlap area 

S = TT^J {x^){y^) of the reaction zone; these represent the 
corresponding mean values of events in this centrality 
class. 

The KLN model involves a couple of parameters that 
need to be adjusted to obtain the correct final charged 
hadron multiplicity dNch/drj in central Au-(-Au collisions 
at RHIC. In [32] this adjustment was performed for ideal 
fluid dynamics (which conserves entropy) coupled to a 
hadron cascade. The model then correctly predicts the 
charged hadron multiplicities at all other collision cen- 
tralities. In our viscous hydrodynamic model, viscous 
heating produces additional entropy, leading to some- 
what larger final multiplicities. We thus perform an 
iterative renormalization of the initial entropy density 
profile until the measured charged hadron multiplicity 
in the 0—5% most central 200^ GeV Au-|-Au collisions 
at RHIC is once again reproduced. The lower panel of 
Fig. 1 shows that, after this renormalization, the model 
again correctly describes the measured [34] centrality de- 
pendence of charged hadron production. The centrality 
dependence of viscous entropy production (which is rela- 
tively larger in peripheral than in central collisions [24]) 
is (at least at RHIC energies) sufficiently weak to not 
destroy the agreement of the model with experimental 
observations. 

The ability of the MC-KLN model to describe the cen- 
trality dependence of charged hadron production with- 
out additional parameters is the main reason for choos- 
ing it over the MC-Glauber model as our basis for ex- 
trapolation from RHIC to LHC energies. It was re- 
cently shown [35] that this centrality dependence is ro- 
bust against running coupling corrections [36-38] in the 
Balitsky-Kovchegov evolution (on which the KLN model 
is based) which were found to hardly affect its shape. 
They do, however, modify the collision energy depen- 
dence of particle production, with the LHC Pb-t-Pb data 



http: / /www. aiu. ac . jp/~ynara/inckln/. 
^ Note that about 10% larger overlap areas are obtained when 
using the entropy density as weight [10, 18], whereas for all but 
the most central collisions the eccentricities of the energy and 
entropy densities are nearly identical [33]. 
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FIG. 1; Centrality dependence of the charged hadron multi- 
plicity per unit pseudo-rapidity, dNch/dri/{NpB,rt/2) as a func- 
tion of iVpart, in 200y4GeV Au+Au collisions at RHIC (bot- 
tom panel) and in (2.76 — 5.5) ATeV Pb-|-Pb collisions at the 
LHC (top panel). Experimental data are from the PHOBOS 
Collaboration [34] for Au-|-Au collisions at -/s = 200 y4 GeV, 
and from the ALICE Collaboration [39] for Pb-|-Pb collisions 
at ^ = 2.76ATeV. The lines are from the MC-KLN model 
(see text). For Au-|-Au at RHIC the MC-KLN model was nor- 
malized to the measured multiplicity in the 0-5% centrality 
bin; at the LHC, the lines bounding the shaded region were 
normalized to diVch/d?7 = 1548 and 1972 (or dAf^h/dy = 1800 
and 2280) at 0-5% centrality, respectively. 



being better described if running coupling corrections are 
included [39] . Our version of the MC-KLN model does 
not include running coupling corrections,^ and we must 
normalize the initial entropy density profile for Pb-|-Pb 
collisions at the LHC separately from Au+Au collisions 
at RHIC. Without such an independent renormalization, 
we overpredict the measured charged multiplicity from 
central 2.76 ATeV Pb+Pb collisions [39, 40] by about 
10%. 



After renormalization we obtain the solid lines bound- 
ing the shaded region in the upper panel of Fig. 1, 
with the lower (upper) bound corresponding to Pb-|-Pb 
collisions at 2.76 (5.5) AGcV, respectively. The data 
in that panel are from the ALICE Collaboration for 
Pb+Pb at 2.76 A GeV [39, 40]. (For 5.5AGeV Pb+Pb 
collisions we assumed dNch/dy = 2280 (corresponding to 
dNch/drj~ 1972), based on an extrapolation of Fig. 3 in 
Ref. [40].) One sees that, even without running cou- 
pling corrections, but including viscous entropy produc- 
tion, the MC-KLN model does a good job in describing 
the measured centrality dependence of charged hadron 
production in Pb+Pb collisions at the LHC. This gives 
hope that the successful description of the centrality de- 
pendence of hadron spectra and elliptic flow at RHIC 
energies (see below and [10]) translates into a reliable 
prediction of the corresponding centrality dependences 
in Pb+Pb collisions at the LHC. 

Figures 2 and 3 establish our baseline for the ex- 
trapolation to LHC energies. In Fig. 2 we show our 
purely hydrodynamic fit (obtained with parameters tq, 
77/s, Tchcm, and Tdoc set as described above^) of the 
hadron spectra measured in 200 A GeV Au+Au collisions 
at RHIC. Fig. 2a shows the mid-rapidity transverse mo- 
mentum spectra per unit pseudo-rapidity for unidentified 
charged hadrons from the STAR [41] and PHENIX [42] 
experiments compared with the hydrodynamical model. 
Figs. 2b, c show a similar comparison for the pr-spectra 
per unit rapidity of identified pions and protons from 
STAR [43, 44] and PHENIX [45]. In the experimental 
spectra, protons from weak decays were removed; STAR 
quotes a large systematic error associated with this feed- 
down correction, and within that large error band the 
two data sets agree with each other, even if the central 
values of the STAR proton data appear to be up to 50% 
higher than PHENIX data. Our results agree well with 
the STAR protons for px > 0.6GeV/c but overpredict 
the PHENIX protons by up to a factor 2. 

Figure 3 shows the hydrodynamically calculated dif- 
ferential elliptic flow for unidentified charged hadrons 
in comparison with STAR ?;2{4}(pt) data [46, 47], for 
four centrality classes ranging from semi-central to mid- 
peripheral collisions (10—50% centrality). With 77/s = 
0.20, viscous hydrodynamics gives an excellent descrip- 
tion of the STAR W2{4} data, even up to 3GeV/c in 
transverse momentum (i.e. beyond the pr range where 
the hydrodynamic description is expected to begin to 



* The rcBK code in [35] includes running coupling corrections but 
it has not been renormalized to take into account viscous entropy 
production at RHIC energies. 



Note that our value tq =0.6fm/c is 45% smaller than the value 
of 1.05 fm/c used for 77/5 = 0.2 in the VISHNU simulations in [10]. 
The earlier evolution of hydrodynamic transverse flow arising 
from this smaller tq value compensates for the lack of a highly 
dissipative hadronic phase in the purely hydrodynamic approach. 
Hadronic dissipation leads to a significant broadening in partic- 
ular of the proton py-spectra during the hadronic stage which 
(given the constraints from the elliptic flow data which pro- 
hibit us from simply lowering T^cc) viscous hydrodynamics with 
temperature-independent = 0.2 cannot replicate. 
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FIG. 2: pT-spectra of all charged hadrons (a), positive pions 
(b) and protons (c) for 200 A GeV Au+Au collisions of dif- 
ferent centralities as indicated. The symbols show data from 
the STAR ([41, 43, 44], x) and PHENIX ([42, 45], +) exper- 
iments. The lines are results from the viscous hydrodynamic 
model for constant = 0.20 and MC-KLN initial conditions 
(see text for other model parameters). 



break down, due to the increasing influence of hard pro- 
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FIG. 3: (Color online) Differential elliptic flow V2{pt) for 
charged hadrons from 200 A GeV Au-|-Au collisions of differ- 
ent centralities as indicated. Open symbols are experimental 
data from the STAR experiment for i;2{4}(pt) [46, 47], lines 
with filled symbols of the same shape are the corresponding 
hydrodynamic fits with the same model as in Fig. 1. For the 
40—50% centrality bin data and theory are vertically offset 
by 0.1 for better visibility. 



duction processes and large uncertainties in the viscous 
correction 5f to the local phase-space distribution at ki- 
netic freeze-out [48]). Looking more carefully one sees 
that our model slightly overestimates the elliptic flow at 
low < 1 GeV while underestimating it in the high-p^ 
region, pt > 2 GeV. 

In [26] wc noted a tension in trying to simultaneously 
fit within a purely viscous hydrodynamic approach the 
proton pT-spectra and the charged hadron differential el- 
liptic vf'{2}{pT) when using EOS s95p-PCE. Even a tem- 
perature dependent rj/ s{T) that has a larger shear viscos- 
ity in the hadronic phase could not resolve this tension: 
in [48] two of us found that the RHIC Au-|-Au hadron 
spectra are insensitive to a temperature-dependent in- 
crease of the shear viscosity in the hadron gas phase, 
as was previously seen in [49]. Figs. 2 and 3 demon- 
strate that this problem is largely resolved when using 
the W2{4}(pt) data (Fig. 3) instead of f2{2} (see Fig. 8 
further below) : We obtain an excellent description of the 
differential elliptic flow, together with an acceptable de- 
scription (within large experimental uncertainties) of the 
Pt spectra. 

Overall, the viscous fluid dynamic description of the 
hadron spectra and charged hadron elliptic flow V2{pt) 
shown here is of similar quality as the hybrid model de- 
scription with VISHNU presented in [10]. Since purely hy- 
drodynamic simulations are numerically much less costly 
than calculations with VISHNU, we will now use them to 
generate a broad range of predictions for soft hadron pro- 
duction in Pb+Pb collisions at the LHC. 
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III. PREDICTIONS FOR PB+PB COLLISIONS 
AT THE LHC 

As discussed above, the extrapolation from RHIC to 
LHC is done keeping tq, Tchem, ^doc and 77/ s fixed. When 
comparing the resulting viscous hydrodynamic predic- 
tions with experimental data from the recently started 
LHC heavy-ion collision program, wc will search for in- 
dications from experiment that would motivate chang- 
ing these parameters. First results for p^-spectra [50] as 
well as both the pT-differential and pr-integrated elliptic 
flow of unidentified charged hadrons [1] have already been 
published and will be compared with the theoretical pre- 
dictions below. Additional experimental information on 
spectra and elliptic flow of identified hadrons will become 
available soon; the relevant hydrodynamic predictions are 
presented in this section. 




Pj (GeV) 

FIG. 4: (Color online) pT-spectra of all charged hadrons, pos- 
itive pions, and protons for minimum bias 200 A GeV Au+Au 
(thin red lines and data points) and (2.76-5.5) A TeV Pb-|-Pb 
collisions (black lines with shaded area). The RHIC data are 
from the PHENIX experiment [45]. The shaded bands for 
the LHC predictions are Umited at the bottom (top) by lines 
for Vs = 2.76 (5.5) ATeV, corresponding to diVdi/dy = 1800 
(2280) (diVch/d?7 = 1548 (1972)). The calculations assume the 
same constant r;/s = 0.2 at all shown collision energies. 

In Fig. 4 we show the transverse momentum spectra 
for all charged hadrons, as well as for identified pions 
and protons, for minimum bias collisions of Au+Au at 
RHIC and Pb-hPb at the LHC.*' For RHIC we compare 
with data from the PHENIX CoUaboration [45]. The 
upper and lower bounds of the shaded areas are predic- 
tions for minimum bias Pb-(-Pb collisions at collision en- 



ergies of 5.5 and 2.76 TeV per nucleon pair, respectively. 
The LHC spectra are visibly flatter than at RHIC en- 
ergies, reflecting stronger radial fiow. For central col- 
lisions (0—5% centrality), the fireball lifetime increases 
from Au-^Au at RHIC to Pb-hPb at LHC by about 19% 
and 24%, respectively, for 2.76 and 5.5 A TeV collision en- 
ergy; for peripheral collisions at 70—80% centrality, the 
corresponding lifetime increases are even larger (34% and 
41%, respectively). The average radial flow velocity in- 
creases in central collisions (0—5% centrality) by 5 and 
7%, respectively, and in peripheral collisions (70—80% 
centrality) by 9 and 11%. 
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FIG. 5: (Color online) pT-integrated elliptic flow of charged 
hadrons for 200 A GeV Au-|-Au collisions at RHIC (open sym- 
bols are STAR data [46, 47], the lower red line is the result 
from viscous hydrodynamics) and for 2.76 ATeV Pb-|-Pb col- 
lisions at the LHC (filled symbols are ALICE data [1], the 
upper magenta line shows the viscous hydrodynamic predic- 
tion). In both experiment and theory the differential ellip- 
tic flow 1)2 (pt) (see Figs. 3 and 7) was integrated over the 
range 0.15 GeV/c<pT < 2 GeV/c for Au-hAu at RHIC and 
over pT > 0.2 GeV/c for Pb+Pb at the LHC. 

Figure 5 shows the integrated charged hadron elliptic 
flow V2 as a function of collision centrality for Au+Au 
collisions at RHIC and Pb+Pb collisions at the LHC. At 
RHIC energy, our results (lower red line) overestimates 
the STAR W2{4} data by about 11% in mid-central colli- 
sions, but agrees nicely with 1)2 {EP} except for the most 
peripheral collisions.^ At first sight the overprediction 
of the PT-integrated W2{4} at RHIC is surprising, given 



To simulate minimum bias collisions, we compute the spectra 
for the centrality classes shown in Figs. 2(b) and 6 and average 
them. Any additional observables, such as the minimum bias 
elliptic flow in Fig. 8 below, are calculated from those minimum 
bias spectra. 



In very peripheral collisions, the fireball lifetime decreases dra- 
matically, cutting short the buildup of anisotropic hydrodynamic 
flow and thereby prohibiting V2 from saturating. In addition, 
viscous effects are stronger in the small fireballs created in pe- 
ripheral collisions than in the larger central collision fireballs. 
Both effects together cause the theoretical V2 values to decrease 
sharply at large collision centralities, in apparent confiict with 
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FIG. 6: pT-spectra per unit pseudorapidity for charged hadrons (a) and per unit rapidity for pious (b) and protons (c) for 
Pb+Pb collisions at the LHC. The lower and upper end of the shaded bands represent viscous hydrodynamic predictions for 
^/s = 2.76 and 5.5 A TeV (corresponding to dNch/drj = 1548 and 1972, or dNch/dy = 1800 and 2280), respectively. Experimental 
data in panel (a) are from the ALICE experiment [50]. 



the excellent description of the differential elliptic flow 
D2{4}(pt) shown in Fig. 3. The apparent paradox is re- 
solved by observing that the hydrodnamically computed 
charged hadron p^-spectra shown in Fig. 2 are somewhat 
harder than measured, thereby giving too much weight in 
the pT-integral to the range 0.75 < pr < 2GeV/c where 
W2{4}(pt) is large. ^ 

At LHC energy (y^^ 2.76 A TeV) our integrated V2 



the experimental data. The experimental V2{2} and i'2{EP} 
measurements are, however, contaminated by non-flow effects, 
in particular in very peripheral collisions. Once non-flow effects 
are corrected for [51], the experimental V2 values decrease at 
large collision centralities much in the same way as predicted by 
hydrodynamics. 

The agreement with the ■U2{EP} data is fortuitous and should, in 
fact, not happen since the measured i'2{EP} includes a positive 
contribution from event-by-event V2 fluctuations [52] while our 
hydrodynamic calculation yields the average elliptic flow (-02) 
which is smaller. 



lies between W2{2} and 'y2{4} values measured by the 
ALICE Collaboration [1]. Again, we overpredict the 
PT-integrated V2{4} by about 10—15%. We note that 
from RHIC to LHC the hydrodynamically computed in- 
tegrated V2 in mid-central to mid-peripheral collisions in- 
creases by about 30%, in agreement with the experimen- 
tal observations. This is due to reduced viscous suppres- 
sion effects in the larger and denser fireballs created at 
the LHC and a longer fireball lifetime which allows the 
momentum flow anisotropy to approach saturation more 
closely than at lower energies [23, 53]. In very peripheral 
collisions, even at LHC energies such a saturation of V2 
does not happen; this is the reason why in Fig. 5 the inte- 
grated V2 is seen to decrease at large collision centralities, 
both at RHIC and LHC. 

In Fig. 6 we present hadron transverse momentum 
spectra for Pb4-Pb collisions at LHC energies, for a range 
of collision centralities. In panel (a) we compare the 
hydrodynamic predictions with first data from the AL- 
ICE experiment [50]. Overall, the theoretical descrip- 
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FIG. 7: (Color online) Differential elliptic flow V2{pt) for unidentifled charged hadrons (a) and identifled pions (b) and protons 
(c), for Pb+Pb collisions of different centralities at the LHC. Experimental data for charged hadron V2{pt), denoted by solid 
symbols, are from the ALICE experiment [1]; they should be compared with theoretical lines carrying open symbols of the 
same shape and color. The shaded bands show the variation of the hydrodynamic predictions with collision energy between 
-^5 = 2.76 and 5.5 TeV (corresponding to dNch/dy= 1800 and 2280, respectively). The lines corresponding to the lower collision 
energy (-^5 = 2.76 TeV) deflne the lower end of the shaded regions at pr = S GeV/c. 



tion of these experimental data is of similar quality as 
for Au+Au collisions at RHIC (see Fig. 2). In the most 
central collisions, the hydrodynamical model describes 
the charged hadron spectrum somewhat better than at 
RHIC, whereas in the very peripheral collisions the hy- 
drodynamic spectra are too flat, presumably due to large 
viscous shear pressure effects. Future comparison with 
the measured spectra at other collision centralities and 
for identified hadrons, shown here in panels (b) and (c) 
as predictions, should shed further light on the origin of 
the discrepancy in peripheral collisions. 

Figure 7a shows a comparison of the hydrodynami- 
cally generated differential vf^{pT) for charged hadrons 
with the ALICE W2{4} data [1], for four different colli- 
sion centralities. For the most peripheral of these, we 
also show the measured 'y2{2} for comparison. The hy- 
drodynamic predictions agree nicely with the data at low 
Pt < IGeV/c, but overshoot the experimental values 
by 10—20% at larger pr, especially in the more periph- 
eral bins. In the 40—50% centrality bin, the theoretical 
prediction happens to agree nicely with W2{2}(pt) even 
though the latter should be shifted upward by flow fluc- 
tuations that arc not included in the theoretical calcula- 
tion. We note that the theoretical overshoot is less severe 
in the VISHNU hybrid model (see Fig. 3 in [20]) than in 
the purely hydrodynamic simulations shown here. This 



suggests that the excess of V2{pt) over the measured val- 
ues at Pt > 1 GeV/c in Fig. 7a may be caused by an 
inadequate description of the late hadronic stage and its 
freeze-out. 

We can summarize Figs. 2a, 3, 5, 6, and 7a by observ- 
ing that the hydrodynamic model overpredicts the px- 
integrated charged hadron V2 by 10—15% at both RHIC 
and LHC, but for different reasons: at RHIC the differ- 
ential elliptic Sow W2,ch(PT) is correctly reproduced while 
the inverse slope of the theoretical pT-spectra is slightly 
too large, while the LHC pr-spectra are described a bit 
better (at least in the most central collisions where pub- 
lished data are available) but the slope of W2,ch(PT) at the 
LHC is slightly overpredicted. 

Panels (b) and (c) of Fig. 7 give predictions for the dif- 
ferential V2(j>t) of identifled pions and protons. Please 
note the different shape of the proton V2(j>t) from that of 
the pions at low pT- radial flow pushes the proton ellip- 
tic flow to larger values of px- Comparing the curves for 
y/s = 2.76 and 5.5 ^ TeV, we see that this "radial push" of 
the proton V2 increases with collision energy, so for higher 
the rise of V2{pt) is shifted to larger transverse mo- 
menta, while at flxed px < 1.5GeV/c the proton elliptic 
flow decreases with increasing collision energy. This hap- 
pens only for heavy hadrons but not for the much lighter 
pions (see panel (b)). 
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FIG. 8: (Color online) Differential elliptic flow V2{pt) for 
all charged hadrons (a) and identified pions and protons (b), 
for minimum bias 200 A GeV Au+Au collisions at RHIC and 
(2.76-5.5) A TeV Pb+Pb collisions at the LHC. Experimental 
data for W2{2} from Au+Au collisions at RHIC are from the 
STAR experiment [54]. Solid lines are viscous hydrodynamic 
results for 200 A GeV Au+Au collisions with the same hydro- 
dynamic parameters as in Figs. 1-4; note their disagreement 
with the 1^2 {2} data shown here (in contrast to their excel- 
lent agreement with «2{4} data shown in Fig. 3). The shaded 
bands are LHC predictions and show the variation of the the- 
oretical expectations for Pb+Pb collisions at collision ener- 
gies ranging from -^5 = 2.76 to 5.5 A TeV (corresponding to 
dNch/dy = 1800 and 2280, respectively). As in Fig. 7, the lines 
defining the lower end of the shaded region at pT = 3GeV/c 
correspond to the lower LHC energy = 2.76 A TeV. 



trality bins: the heavier the hadron, the stronger a push 
of V2 towards higher pT is observed. At sufficiently large 
Pt, V2{pt) is larger at LHC than at RHIC for all particle 
species, but at low pt this holds only for pions whereas 
all heavier hadrons show a decrease of V2{pt) from RHIC 
to LHC at fixed pr- As the hadron rest mass grows, 
the crossing point where the decrease of V2 at fixed pr 
with rising collision energy turns into an increase shifts 
to larger pT values. In view of Fig. 9, the experimen- 
tal observation [1] that for unidentified charged hadrons 
v'^{pt) hardly changes at all from RHIC to LHC appears 
accidental:^ The increase of V2{pt) at fixed pT for pions 
is balanced by a corresponding decrease for all heavier 
hadrons leaving, as it happens, no visible net effect once 
all charged hadrons are lumped together. 

In Refs. [10] it was argued that a robust method for 
extracting the QGP shear viscosity is to fit the collision 
centrality dependence of the eccentricity-scaled charged 
hadron elliptic flow / ewith a viscous hydrodynamic + 
hadron cascade hybrid code. In that study it was found 
that, at fixed collision energy,^° plotting vf^/e against 
the charged hadron multiplicity density per unit overlap 
area, {l/S){dNch/dy), yields "universal" curves that de- 
pend only on the QGP shear viscosity but not on the 
model for the initial energy density distribution (in par- 
ticular its eccentricity). In Fig. 10 we show such a plot 
for 200 A GeV Au+Au collisions at RHIC together with 
Pb+Pb collisions at two LHC energies. The four pan- 
els show this scaling in terms of distributions in pseudo- 
rapidity {rj, left column) or rapidity (y, right column), 
and also compare it for our default choice of using the 
initial energy density as weight for the calculation of 
the average eccentricity s and overlap area 5* (top row) 
with what one obtains by evaluating these quantities 
with the initial entropy density instead (as is done in 
Refs. [10, 18]) (bottom row). We sec that, independent of 
these choices of representation, the universality of vf^/e 
vs. {1/ S){dNch/dy) or {l/S){dNch/dri) does not carry 
over to different collision energies (at least not for the 
purely hydrodynamic simulations studied in the present 
work): At the same multiplicity density {1/ S)(dNch/dy) 
or {\/ S){dNf±_l drf), more peripheral higher energy col- 
lisions produce less elliptic flow per initial eccentric- 
ity than more central lower energy collisions. At fixed 
i{/s = 0.2, the difference between 200 A GeV Au+Au and 



In Figs. 8 and 9 we pursue this theme further, by di- 
rectly comparing the differential elliptic fiows at RHIC 
and LHC energies. In Fig. 8 we show results for minimum 
bias collisions; the RHIC predictions are compared with 
available data from STAR [54]. We see that at low px, 
the elliptic flow for unidentiflcd charged hadrons (which 
are strongly pion dominated) and for identified pions in- 
creases from RHIC to LHC whereas the opposite is true 
for protons. At higher [pT > 1.5 GcV/c), on the other 
hand, V2{pt) increases for both pions and protons as we 
increase the collision energy. Fig. 9 shows this for a few 
more hadron species, for the 10—20% and 40—50% cen- 



^ Contrary to the claim made in [12], the observation that the 
ratio between ^^'^(pt) measured at LHC and at RHIC is approx- 
imately independent of pT cannot be directly used to conclude 
that (»7/.s)qgp does not change from RHIC to LHC. If that ar- 
gument were correct, this ratio should be independent of pj- not 
only for the sum of all charged hadrons, but also for each iden- 
tified hadron species separately. Our hydrodynamic calculations 
show that the latter does not hold even if remains unchanged 
from RHIC to LHC. 

" Wc recently checked that this multiplicity scaling carries over 
to other collision systems such as Cu+Cu at the same collision 
energy [55]. 
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FIG. 9: (Color online) Comparison of the differential elliptic flow V2{pt) for 200^ GeV Au+Au collisions at RHIC (dashed 
lines) and 2.76 A TeV Pb+Pb collisions at the LHC (solid lines), at 10%-20% (a,b) and 40%-50% (c,d) centrality, for a variety 
of different hadron species. Note the slightly negative elliptic flow for the heavy hyperons at low pr- 



2.76 A TeV Pb+Pb collisions (red circles vs. green up- 
ward triangles in Fig. 10) is as large as the difference 
between 77/5 = 0.16 and 77/5 = 0. 20 for Au+Au collisions 
at fixed -^5 = 200 A GeV (red circles vs. black squares). 

We note that the tendency in Fig. 10 of higher energy 
collisions producing less vf^/e at fixed {1/ S){dNch/dy) 
than lower energy ones contradicts the opposite ten- 
dency observed by Hirano et al. in Fig. 3 of Ref. [18] 
where an ideal hydro + hadron cascade hybrid code was 
employed. The authors of [18] presented strong ar- 
guments that their observation of larger v^/e at fixed 
{l/S){dNch/dri) in higher energy collisions is not related 
to their use of a hadron cascade for describing the late 
hadronic stage. Our opposite finding, on the other hand, 
is supported by the earlier purely hydrodynamic scaling 



The careful reader will notice that for 200 A GeV Au-|-Au col- 
lisions, our maximal values for (1/ S){dN^i^/dr]) in Fig. 10b are 
significantly larger than those shown in Fig. 3 of Ref. [18] . This is 
due to a lower normalization of the initial entropy density in [18], 
corresponding to dN^^/drj ~ 600 instead of our dN^^/dr) ~ 700 
in central Au+Au collisions (T. Hirano, private communication). 



studies presented in the last two works of [24] whose au- 
thors came to the same conclusion as we do here. At 
present this discrepancy remains unresolved; we suspect, 
however, that the origin of the difference between our 
work and that of Hirano et al. could be in their use of 
a more realistic (3+l)-dimensional hydrodynamic evo- 
lution [56], although in the earlier ideal fluid hydrody- 
namical studies at the full RHIC energy, the differences 
between boost invariant and non-boost invariant results 
were small [3, 57]. Possible consequences of the violation 
of boost-invariance in RHIC and LHC heavy-ion colli- 
sions are presently being studied [55]. 

Before moving on, let us comment on the different 
shape at the high-multiplicity end of the curves shown in 
Fig. 10 when using entropy instead of energy density as 
the weight for calculating the initial eccentricity e over- 
lap area S: It is caused by the different centrality de- 
pendence of the energy and entropy density weighted ec- 
centricities in near-central collisions observed in Ref. [33] 
whose authors showed that in the most central collisions 
(where e is dominated by event-by-event shape fluctu- 
ations) the entropy- weighted participant eccentricity de- 
creases faster with decreasing impact parameter than the 
energy-weighted one. 
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FIG. 10: (Color online) Eccentricity-scaled elliptic flow v-^js as function of the charged hadron multiplicity density per unit 
overlap area S from viscous hydrodynamic calculations at -^ = 0.2, 2.76, and 5.5ylTeV (corresponding to dA'ch/dy = 814 
(Au+Au), 1800 and 2280 (Pb+Pb), respectively). Each line corresponds to one collision system at fixed collision energy but 
different coUision centralities (from right to left, the symbols correspond to 0-5%, 5-10%, 10-15%, 15-20%, 20-30%, 30-40%, 
40-50%, 50-60%, 60-70% and 70-80% centrality). The four panels show ^;2(r;)/£ vs. (l/S')(diVch/d??) (where r] is pseudorapidity) 
(a,c) and V2{jj)/e vs. {l/S){dNch/dy) (where y is rapidity) (b,d), with e and S evaluated with the participant-plane averaged 
energy density e{r±,To) as weight function (a,b) (default option, see Sec. II) or (for comparison with other work) with the 
corresponding entropy density s{r±,To) as weight (c,d). The RHIC curves for 77/5 = 0.16 (black squares) illustrate the effect of 
changing the value of the specific shear viscosity by 0.04. The LHC calculations are done for 77/5 = 0.20 as obtained from the 
hydrodynamic fit to RHIC data. 



IV. TEMPERATURE DEPENDENT r]/s{T) 



explored in this section: 



Shear viscosity is known to suppress the buildup of 
eUiptic flow. Naively, the systematic overprediction of 
V2{4:}{pt) m Pb+Pb collisions at the LHC seen in Fig. 7a, 
together with the excellent description of the same quan- 
tity in Au+Au collisions at RHIC seen in Fig. 3, thus 
suggests that the fireball liquid formed in LHC collisions 
might be slightly more viscous (i.e. possess larger aver- 
age 77/5) than at RHIC energies [20, 49]. In this section 
we present some results using a temperature dependent 
specific shear viscosity, (-)(T), that were motivated by 
such considerations. 

Figure 11 illustrates the following three trial functions 



^) = 0.2 + 0.3 

cm 

= 0.2 + 0.4.^^"^'='^""^' 

s / 2 



rp2 

chcn 



'A = 0.2 + 0.3 



s/ 3 



chcn 



(1) 

(2) 
(3) 



Here Tchom = 165 McV is the chemical decoupling tem- 
perature and stands for the "transition temperature" at 
which the hadronization of quarks and gluons is com- 
plete. 

In principle, the value of 77/s should exhibit a minimum 
near Tchem and increase again in the hadronic phase be- 
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FIG. 11: (Color online) Three temperature dependent 
parametrizations {ri/s){T) studied in this section. In all cases, 
77/s = 0.2 for T < Tchem = 165 MeV. 



low Tchcm [58-60]. The authors of [49] pointed out, how- 
ever, that at the full LHC collision energy of S.SATeV 
the behavior oirj/s at temperatures below Tchcm has very 
little effect on the final hadron spectra and their elliptic 
flow. At 2.76 A TeV the effect on elliptic flow was mod- 
erate, and negligible on the spectra. Here we will con- 
centrate on qualitative aspects of effects arising from a 
temperature dependent growths of 77/s in the high tem- 
perature region that can be explored at LHC energies 
but is beyond the reach of RHIC, and continue to set 
rj/s = 0.2 at T < Tchem for simplicity. 

As pointed out in [49], the spectra and elliptic flow in 
Au-|-Au collisions at RHIC energies are most sensitive to 
the average value of 77/s in the temperature region be- 
low 220-230 MeV. We have checked that ahering 7;/s at 
higher temperatures as shown in Fig. 11 has little influ- 
ence on the results at RHIC energies shown in Sec. II. 

Figure 12 illustrates the influence of a linear tem- 
perature dependence of 77/s as in Eq. (1) on the cen- 
trality dependence of charged hadron production. The 
solid black line is the same as shown in the upper 
part of Fig. 1 where it forms the lower bound of the 
shaded region; it corresponds to constant ?7/s = 0.2 and 
Navier-Stokcs initial conditions for the shear stress ten- 
sor, tt'^'^ = at tq =0.6fm/c. The dashed and dash- 
dotted lines in Fig. 12 use {ri/s)i{T) with cither Navicr- 
Stokes (dashed) or zero (dash-dotted) initial conditions 
for TT^''. These last two lines were normalized to the 
ALICE point for the 0-5% most central Pb-FPb colli- 
sions (dA^ch/rf?7= 1584±80 [40]), whereas the black line 
was normalized to our best guess before the ALICE 
data became available {dNch/dy = 1800, corresponding 
to (iiVch/'i77 = 1548). The centrality dependence is then 
controlled by the predictions from the MC-KLN model, 
modified by viscous entropy production during the hy- 
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FIG. 12: (Color online) Final charged multiplicity per 
wounded nucleon pair as a function of number of participant 
nucleons in Pb-|-Pb collisions at v^ = 2.76ATeV, for differ- 
ent functional forms of {rj/s){T) and initial conditions for the 
shear stress tensor vr''" (see text). 

drodynamic evolution. 

We see that even a relatively modest temperature de- 
pendent increase of 77/s in the QGP phase leads to a 
significantly stronger non-linearity in the dependence of 
charged particle production on the number of wounded 
nucleons. The reason is that an increase of 7;/s with 
temperature leads to more viscous heating in central col- 
lisions (which probe higher initial temperatures and such 
larger effective shear viscosities) than in peripheral ones 
(whose initial temperatures are lower). Since the entropy 
production rate is given by 



this effect is stronger for Navier-Stokes initial conditions 
(where tt^'^ is proportional to the velocity shear tensor 
CT^'' which at early times diverges like 1 /r) than for zero 
initial shear stress (where tt'"^ starts from zero and ap- 
proaches its Navier-Stokes value 2ria'^'^ only after several 
relaxation times when, due to its l/r decay, it has 
already decreased to much smaller values). 

If one were to postulate the validity of the MC-KLN 
model as the correct description of the initial particle pro- 
duction, the ALICE data shown in Fig. 12 would exclude 



For reference we list the fractions of the finally measured entropy 
in the most central and most peripheral centrality bins shown in 
Fig. 12 that are generated by viscous heating during the hydrody- 
namic expansion: Constant rj/s = 0.2: AS/Sflnal = 26% (0—5%) 
and 33% (70-80%); (r;/s)i(r) with ir^'" = 0: AS/Sflnal = 25% 
(0-5%) and 15% (70-80%); {r)/s)i{T) with tt^" = 2r)o-'"^ : 
AS/Sflnai = 60% (0-5%) and 49% (70-80%). 
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a temperature dependence of 77/s as given in Eqs. (1) and 
(2) for Naver-Stokes initial conditions. While we are not 
prepared to make such a statement on the basis of Fig. 12 
alone, we believe that it is important to point out this 
relatively strong sensitivity of the centrality dependence 
of dNch/dr] to the transport properties of the expanding 
fireball medium and to emphasize the constraints it thus 
places on possible models for the QGP shear viscosity. 
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FIG. 13: (Color online) Charged hadron transverse momen- 
tum spectra (top) and differential elliptic flow (bottom) for 
2.76 yl TeV Pb+Pb collisions at 20-30% centrality, for differ- 
ent models for the temperature dependence of 77/s and dif- 
ferent initial conditions for tt'"^ (Navier-Stokes ("NS") or 
("Zero")). The ALICE data in the bottom panel are from 
Ref. [1]. 



We now turn to the discussion of the influence of a 
possible temperature dependence of 77/s on the charged 
hadron p^-spectra and elliptic flow. Figure 13 shows 
LHC predictions for 2.76 A TcV Pb+Pb collisions of 
20—30% centrality. To ensure comparability of the dif- 
ferent cases studied in this figure we simply normal- 



TABLE I: Initial central entropy densities so and tempera- 
tures To for the viscous hydrodynamic simulations of 20—30% 
centrality Pb-|-Pb collisions at the LHC (ys = 2.76 ^ TeV) 
shown in Fig. 13. The different models for the T-dependence 
of 7;/s are defined in Eqs. (1)~(3). "0" stands for n^'^ = at 
To, "NS" stands for Navier-Stokes initialization of the shear 
stress tensor, ttq" — 2r]a'^'' at tq. 



77/s model 


K) 


so (fm-^) 


To (MeV) 


r7/s = 0.2 




NS 


191.6 
172.4 


427.9 
413.9 






NS 


179.6 
119.3 


419.2 
368.7 


{n/s)2{T) 




NS 


179.6 
115.6 


419.2 
365.1 


{n/s)z{T) 




NS 


175.2 
116.6 


416.0 
366.1 



ized the initial entropy density profile such that we al- 
ways obtain dNch/dy = 770, i.e. the same value that we 
had obtained before for constant 77/5 = 0.2 at this cen- 
trality. We first note that for constant 77/5 = 0.2, we 
don't observe any significant difference in the charged 
hadron spectra and elliptic flow between zero and Navier- 
Stokes initialization for tt'"' . Turning to the temperature- 
dependent parametrizations [rj/ s)i{T), wc note that for 
zero initialization of ■k^'^ (solid lines) our results agree 
with those reported in [49]: An increase of 77/s at higher 
QGP temperatures leads to somewhat harder charged 
hadron py-spectra (i.e. somewhat stronger radial fiow, 
caused by the larger tranverse effective pressure gradi- 
ents at early times) and a suppression of the differen- 
tial elliptic flow (due to an increase of the timc-averaged 
effective shear viscosity of the fluid). It is interesting 
to observe the hierarchy of the curves in Fig. 13 corre- 
sponding to the three parametrizations (l)-(3): For the 
PT-spectra, all three T-dependent viscosities lead to al- 
most identical hardening effects on the spectral slope, 
while for the differential elliptic flow v^{pt) the curves 
are ordered not according to the 77/s-values at the initial 
central fireball temperature (see Table I), but according 
to their hierarchy in the 165<T<280McV range. In 
fact, the observed magnitudes of the viscous V2 suppres- 
sion for the three [rj/ s){T) fimctions suggest that, at this 
beam energy and collision centrality, the buildup of ellip- 
tic flow is dominated by the QGP transport properties 
at 200<r<250McV. (At RHIC energies, the transport 
properties for T < 200—220 McV dominate the generation 

of 772 [49].) 

For Navier-Stokes initial conditions (dashed lines in 
Fig. 13), the increase in radial flow caused by an in- 
crease of 77/5 at high temperature is stronger and the 
viscous V2 suppression is weaker than for zero initial ■n^'^ . 
This is caused by the much larger initial shear stress ten- 
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sor components in the NS case, compared to the case 
of ttq'^ = where tt^'^ approaches its (by that time al- 
ready much smaUer) Navier-Stokes hmit only after sev- 
eral relaxation times [24]. The increase oi r]/s with 
temperature generates a steeper initial transverse effec- 
tive pressure gradient (since tt'"^ grows faster than the 
entropy density s when rj/s increases with temperature), 
and this generates stronger radial flow. It also causes a 
larger spatial eccentricity of the initial effective pressure 
profile which (when compared to the case of tTq = 0) gen- 
erates stronger elliptic flow. In fact, we found that for 
earlier starting times tq (where the Navicr-Stokcs values 
for tt'^" are even larger), the quadratic paramctrization 
{vl s)2{T) with NS initial conditions can lead to more el- 
liptic flow than a constant "(]!& = 0.2, in spite of the larger 
mean viscosity of the fluid. 

We conclude from this exercise that a firm determi- 
nation whether or not the ALICE data point towards a 
temperature-dependent growth of 77/5 with increasing T, 
as expected from pcrturbative QCD [61] and (perhaps) 
from lattice QCD [62], is not possible without a bet- 
ter understanding of the initial conditions for the energy 
momentum tensor (in particular the shear stress compo- 
nents) at the beginning of the hydrodynamic evolution. 
Whereas generically larger viscosities cause a suppression 
of the elliptic flow, temperature-dependent viscosities can 
influence the initial effective pressure profile and its ec- 
centricity in a way that counteracts this tendency and, 
for some models such as Navier-Stokes initial conditions, 
can even overcompensate it. 



ration. Our extrapolation from RHIC to LHC energies 
was based on the assumption that the QGP shear vis- 
cosity 77/ s does not change with increasing fireball tem- 
perature and stays fixed at the value ry/s = 0.2 extracted 
from the RHIC data, assuming MC-KLN initial condi- 
tions. The start time tq for the hydrodynamic evolution 
and the freeze-out temperature Tdcc were held fixed, too. 
We found that, using the beam energy scaling implicit 
in the MC-KLN model, such an extrapolation gives a 
good description of the centrality dependence of charged 
hadron production and the charged hadron py-spectra 
in central Pb-Pb collisions, but overpredicts the slope 
of the pT-differential elliptic flow and the value of its 
PT-integrated value by about 10—15% in mid-central to 
mid-peripheral collisions. In the most peripheral colli- 
sions, the predicted charged hadron p^-spectra are too 
flat, and the integrated elliptic flow is too small com- 
pared to the experimental data. A preliminary study of 
possible temperature dependent variations of 77/5 in the 
high-temperature region explored for the first time at the 
LHC remained unconclusivc but pointed to a clear need 
for better theoretical control over the initial conditions 
for the hydrodynamic energy-momentum tensor, in par- 
ticular its shear stress components. The development of 
detailed dynamical models for the pre-thermal evolution 
of the collision fireball and their matching to the viscous 
hydrodynamic stage is a matter of priority for contin- 
ued progress towards quantifying the transport proper- 
ties of the quark-gluon plasma at different temperatures 
and densities. 



V. CONCLUSIONS 

Based on an successful global fit of soft hadron pro- 
duction data in 200AGeV Au-|-Au colhsions at RHIC 
with a pure viscous hydrodynamic model with Cooper- 
Frye freeze-out, presented in Sec. II, we generated hydro- 
dynamic predictions for the pr-spectra and differential 
elliptic flow of unidentifled charged hadrons and identi- 
fied pions and protons for Pb-|-Pb collisions at the LHC. 
Where available, these predictions were compared with 
available experimental data from the ALICE Collabo- 
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